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ABSTRACT 

Presented are the hrst interferometric images of cool starspots on the chromo- 
spherically active giant A Andromedae. These images represent the hrst model- 
independent images of cool starspots on a star other than the Sun to date. The 
interferometric observations, taken with the Michigan Infra-Red Gombiner cou¬ 
pled to the Genter for High Angular Resolution Astronomy Array, span 26 days 
from Aug 17*^, 2008 to Sep 24*^, 2011. The photometric time series acquired 
at Fairborn Observatory spanning Sep 20*^, 2008 to Jan 20*^, 2011 is also pre¬ 
sented. The angular diameter and power law limb-darkening coefficient of this 
star are 2.759 ± 0.050 mas and 0.229 ± 0.111, respectively. Starspot proper¬ 
ties are obtained from both modeled and SQUEEZE reconstructed images. The 
images from 2010 through 2011 show anywhere from one to four starspots. The 
measured properties of identical starspots identihed in both the model and recon¬ 
structed images are within two a error bars in 51% of cases. The cadence in the 
data for the 2010 and 2011 data sets are sufficient to measure a stellar rotation 
period based on apparent starspot motion. This leads to estimates of the rotation 
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period (P 2010 = 60 ± 13 days, P 2011 = 54.0 ± 7.6 days) that are consistent with 
the photometrically determined period of 54.8 days. In addition, the inclination 
and position angle of the rotation axis is compnted for both the 2010 and 2011 
data sets; valnes ('P = 21.5°, i = 78.0°) for each are nearly identical between the 
two years. 

Subject headings: stars:starspots — starsdmaging — stars:rotation — stars:individnal(A^^H 
And) — methodsdnterferometric — methods:photometric — instrnmentationdnterferometers 
— techniqnesdmage processing 


1. INTRODUCTION 


In the decades since starspots were hrst hypothesized by Kron (1947), starspots have 


been stndied, in detail, on scores of other stars. More recently, space missions, snch as Kepler, 


have increased this number to potentially tens of thousands (Basri et ah 2011). One thing 


that is certain is starspots are an ubiquitous phenomenon found on stars ranging in age from 
pre-main sequence to evolved giants ( Strassmeier]2009 and references therein). A motivation 
for studying starspots is a better understanding of stellar interiors, particularly the origins 
of magnetic dynamos. Another motivation is that starspots complicate measurements of 
fundamental stellar properties (i.e. Teg, L, RV, etc.) Besides the astrophysical importance, 
if a spotted star happens to harbor orbiting exoplanets, the increased uncertainties in the 
stellar properties will translate directly to increased uncertainties in the exoplanet properties 
(i.e. mass, radius). With the advent of millimagnitude photometry, meter per second radial 
velocity surveys and direct milliarcsecond radius measurements, starspots as a “second- 
order” effect can no longer be ignored. 

To account for the affects of starspots on particular measurements, the starspots them¬ 
selves must be properly characterized. The bulk of the current understanding of starspot 
properties (i.e. size, temperature, number, location) stems from two indirect observational 
techniques of magnetically active stars: light curve inversion and Doppler imaging. Both 
methods suffer from certain assumptions that, if incorrect, can lead to unexpected starspot 


properties or imaging artifacts such as polar starspots and/or latitudinal starspot belts (Un- 
ruh & Collier Cameron 1997; Berdyugina & Tuominen 1998). The assumptions include 


a priori knowledge of the stellar inclination for accurate starspot latitude measurements, 
precise information of stellar parameters, accurate stellar atmosphere models, and accurate 
atomic and molecular line lists. While the concern that these are artifacts rather than real 


features has been largely addressed (Unruh 1996 Rice 2002), a more direct method for 


imaging starspots would bolster conhdence in the present results. 
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Laying aside the veracity of these techniques, the starspot characteristics these tech¬ 
niques have provided in terms of lifetimes, sizes, affect on stellar luminosity, and location 
is in many cases contrary to the behavior of sunspots. For instance, large starspots are 


known to persist from months to years (Berdyugina 2005). However, typical sunspot life¬ 


times range from days to weeks. The covering factor, or percentage of the visible surface 
covered by spots, is far larger for active stars (10% to 50%) than for the Sun where the 


covering factor never exceeds 0.2% (Cox 2000). In addition, at times where the covering 


factor is largest, the overall luminosity of active stars decreases substantially (AF < 0.6) 
whereas the Sun’s overall luminosity actually increases. Sunspots appear at a latitude of 
~30° symmetric about the equator at the beginning of the solar activity cycle. As the cycle 
progresses, sunspots migrate toward the equator stopping at a latitude of ~8° (Babcock 


1961, and references therein). Starspots have been observed to reside anywhere from low to 


high latitudes or at the poles (Strassmeier]2009, and references therein). 


Direct imaging of starspots is needed in order to confirm this discrepancy and poten¬ 
tially link the behavior of magnetic spots from solar-type stars to more active stars. This 
direct measure of starspot properties can be obtained via long baseline optical/near infra-red 
interferometry (LBI). By combining the light, akin to Young’s double slit experiment, from 
multiple, widely spaced telescopes, sub-milliarcsecond angular resolutions can be achieved. 


Since the hrst near-IR aperture synthesis image of the binary Capella by Baldwin et ah 


(1996), images of binary stars, rapidly rotating stars, and star-|-disk systems are growing 


more commonplace ( 

Tuthill et al. 

2001 

Hummel 

2004 

Kloppenborg et al. 

2010 

Ghe et al. 

2011 

Baron et al. 

2012). Additionally, bright, convection-induced starspots have been im- 


aged on the late type stars Betelgeuse and VX Sagittarii using LBI (Young et ah 2000 


Chiavassa et al. 2010). More recently similar starspots have also been imaged on the sur¬ 


faces of both RS Per and T Per ( Baron et al.||2014 ). At present, the angular resolution of the 
longest baseline interferometer is ~0.4 milliarcsecond (mas) in the H band. This provides 
the highest imaging resolution by any instrument to date. The median angular diameter 
for a survey of A, F and G main sequence stars is 0.991 mas or ~2.5 resolution elements 
(Boyajian et al. 2012). Therefore the technique is currently only viable for giant stars and 


close early-type dwarfs. 


2. OBSERVATIONS AND DATA REDUCTION 

2.1. The Chromospherically Active Giant A Andromedae 

A Andromedae (A And; HD 222107) is a bright (F: 3.872 mag, H\ 1.501 mag) G8 
IV-HI classihed as a RS GVn type variable star in the Third Gatalog of Ghromospherically 
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Active Binaries ( Eker et al.||2008 ). Calder (1935) first discovered the photometric variability 
of A And with a historical peak AV amplitude of ~0.3 magnitudes. Henry et ah (1995) 


conducted a 15 year photometric monitoring campaign finding a periodic variability of 53.95 


± 0.72 days over an 11.1 ± 0.4 year stellar activity cycle. A And was found by Walker 
(1944) to be a single spectroscopic binary with an orbital period of 20.5212 days. Mass 


estimates for A And in Donati et ah (1995, hereafter D95) range from 0.65-0.85 Mq, and 
with considerable uncertainty. From the mass function provided by the spectroscopic orbit, 
D95 infer a mass ratio q = 0.121qqJ. This in turns yields a companion mass of 0.08-0.10 
Mq leading to the conclusion the companion is a low main-sequence dwarf or high mass 
brown dwarf. The high flux contrast between the two components of A And will preclude 


the companion affecting the photometric or interferometric observations. Nordgren et ah 


(1999) measured a limb darkened angular diameter of A And of 2.66 ± 0.08 mas using the 


Naval Prototype Optical interferometer in the optical; the measurements spanned 10 spectral 
channels ranging from 649 to 849 nm. This angular diameter is a factor of ~5x larger than 
the i7-band angular resolution provided by the CHARA interferometer. In short, A And is 
an interferometrically single, large, bright star with significant variability strongly believed 
to arise from cool starspots. 

Doppler imaging of A And has not been possible due to low projected rotational velocity 
(r’sin(z) = 6.5 km/s) insufficiently broadening absorption lines necessary to detect the de¬ 
formations caused by starspots ( Strassmeier||2009 ). On the other hand, light curve inversion 


has been modestly successful in studying starspots on A And. Frasca et ah (2008, hereafter 


F08) created a surface map via light curve inversion of optical photometry coupled with 
spectral line ratios. This map shows two starspots, each covering ~ 8% of the visible surface 
with temperatures ~880 K cooler than the photosphere, separated by 81° in longitude. In 
addition, each starspot is preceded by a bright active region that is of comparable size to 
the starspot. 


2.2. Interferometric Observations 


For those uninitiated with the terms and physics behind optical/near infra-red interfer¬ 
ometry, an excellent review was written by [Me 


ionnier 


(2003). 


A And was observed on 26 nights between Aug 17*^, 2008 and Sep 24*^, 2011. Table [T] 
lists the date of the observations, the baselines utilized, the number of [u,v] points, and the 
calibrators (dehned below) used on each night. The parenthetical number beside a calibrator 
indicates the number of times it was observed during the night. Table [^contains the uniform 
disk angular diameters with error for each calibrator. All observations were conducted using 
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the Center for High Angular Resolution Astronomy (CHARA) array owned and operated by 
Georgia State University. The array is composed of six 1-m telescopes in a non-redundant 
“Y”-shaped conhguration. The baseline lengths range from 34 to 331 m, currently making 


this the longest baseline optical/near-infrared interferometer in the world (ten Brummelaar 


et al. 2005); the longest baselines provide an angular resolution of ~0.4 mas in the R-Band. 


The data were collected using the image-plane Michigan Infra-Red Combiner (MIRC) in the 
R-band; (see Monnier et ah (2004, 2006) for details). A low resolution (R ~ 42) prism splits 


the light into eight spectral channels with absolute wavelength precision of ±0.25% based 
on measurements of l Peg using the orbit of Konacki et al. (2010). 


Interferometric data are collected when the light path difference between each telescope 
pair is well below the coherence length. This is achieved by “delaying” the light from one 
telescope by adding more path length to the light of the other telescope through the use of 
delay lines. When the light paths are equal, the light from each telescope combines as an 
interference fringe. The MIRC combiner is then set to track these fringes while the data 
frames are taken. Afterwards a series of calibration frames are recorded that include a back¬ 
ground and foreground frame, along with images of the light from each beam individually. 
Collection of the data and calibration frames typically does not exceed 30 minutes. The 
total amount of data taken on a particular night is identical to the number of [u,v] points 
listed in Table [U 


The standard MIRC pipeline was used for data reduction (Monnier et al. 2007). The 
frames containing the fringe pattern in each block of data were co-added. These co-added 
frames are corrected for instrumental effects through a background frame subtraction and 
foreground frame normalization. Raw squared visibilities, triple amplitudes, and closure 
phases are extracted using the Fourier transform of these corrected, co-added frames. Pho¬ 
tometric calibration due to differences in the flux amplitude per telescope beam is performed 
via real-time flux estimates derived from choppers that temporally encode the light from each 
telescope (prior to 2010) or through the use of a beam splitter following spatial hltering to 
shunt part of the beam to a CCD that directly monitors the flux levels from each telescope 
(Che et al. 2010, after 2010). The data are then transformed from relative measurements to 
absolute measurements through observations of a calibration star or “calibrator”. A calibra¬ 
tor is a star of known size that is typically on the order or smaller than the array’s resolution 
limit and is within a few degrees to the target on the sky. 


During the 2008 and 2009, only one or two snapshot observations were obtained of A 
And. A snapshot observation is where only a single set of data frames bracketed by calibrator 
observations is taken. While these observations were valuable conhrming the presence of 
starspots from asymmetric closure phases, the starspots in the resultant images that cannot 
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be confirmed as genuine. The discussion of the observational strategy and results of these 
data sets can be found in Appendix 

The 2010 observing run combined observations from the S1-E1-W1-W2 and S2-E2-W1- 
W2 telescope conhgurations. Each telescope has an alphanumeric designation to describe its 
location. The letter refers to the baseline direction (“S” for south, “E” for east, “W” for west) 
while the number indicates location along that baseline, where “1” is exterior to “2”. The 
longest CHARA baseline is between El and Wl. The switch between telescope conhgurations 
occurred approximately at the mid-point of the night. In both conhgurations A And was 
observed continuously with each observation bracketed by an observation of a calibrator. 
This strategy of telescope conhgurations, constant monitoring and bracketed observations 
is designed to provide well calibrated measurements sampling the largest [u,v] coverage 
available. This strategy yields 11 visibilities, 8 closure phases, and 8 triple amplitudes per 
spectral channel. 

In addition to the above strategy, 2010 data sets consist of observations on sequential 
nights that were subsequently combined into a single OIFITS hleQ This strategy provides 
both an increased [u,v] coverage and a sanity check for the imaging methods since images of 
the star should be nearly identical between subsequent nights. The only exception to this 
strategy is for Sept 10*^ and 11^^, as poor weather prevented observations on the 11^^. 

The 2011 observing run benehted from the MIRC upgrade that enabled simultaneous 
6 telescope observations. Again to maximize the [u,v] coverage, A And was observed con¬ 
tinuously over the night for as long as the delay lines would permit (typically 6 hrs for all 
6 telescopes). Again A And measurements were bracketed by measurements of a calibrator 
star. Each individual observation yields 15 visibilities, 20 closure phases, and 20 triple am¬ 
plitudes per spectral channel. Since observations from subsequent nights were not combined, 
the number of [u,v] points obtained is approximately half of that of the 2010 data set, despite 
the addition of two telescopes. 


2.2.1. Interferometric Measurement Errors 

Standard errors are propagated for the measured quantities. However, because of sys¬ 
tematic effects two types of error are applied to the calibrated squared visibilities and triple 
amplitudes to appropriately account for systematics. Additive errors are necessary for two 
different behaviors in the data. The hrst is when the calibrated squared visibility and/or 


OIFITS is the standard file format for optical/near-IR interferometric measurements (Pauls et al. 


20051. 
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triple amplitude falls below zero. As this is non-physical, a constant is added to enlarge the 
error to include zero. The second is when the squared visibilities and triple amplitudes do not 
monotonically increase or decrease as a function of wavelength across the eight spectral chan¬ 
nels. The errors are then enlarged by a multiplicative constant to account for any abnormal 
structures (e.g. step-functions) found in the data across the 8 spectral channels. Typical 
additive errors for the squared visibility and triple amplitude are 2 x 10“^ and 1 x 10“®, 
respectively. The multiplicative errors improved after 2010 due to better photometric cali¬ 
bration provided by the photometric channels. The typical multiplicative errors in squared 
visibility and triple amplitude are 15% (10%) and 20% (15%), respectively prior to 2010 
(2011). 


Typically 1° is added to the closure phase errors, as suggested by Zhao et ah (2011). In 
addition, two additional closure phase errors are incorporated to avoid poor model fits due 
to calibration systematics. These new errors are important in the low signal-to-noise (S/N) 
regime near to visibility null crossings. As correlated camera readout noise dominates the 
closure phase measurements at low S/N, minimum closure phase errors are applied when the 
S/N in the triple amplitude signal is < 1. Finite time-averaging and spectral bandpass effects 
are accounted for by an error term proportional to ACPa across each spectral channel. CPa 
corresponds to the closure phase as a function of wavelength. These two errors are applied to 
the closure phase noise via the following equation acp > MAX{{30°/{S/N^^^^p),0.2ACPx)), 
where S/ATSamp is the signal-to-noise in the triple amplitude measurement. 


2.2.2. Identification and Removal of Poor Data 

In a few data sets in 2010 and 2011, the visibility in certain blocks of data on certain 
baseline pairs are not consistent with a simple Bessel function. This inconsistency is not 
due to starspots as these features will only change the amplitude of the visibility lobes while 
keeping the form of the function intact. The reason for these discrepant data blocks is most 
likely due to poor calibration on that baseline. These discrepancies are found almost exclu¬ 
sively on the short baseline pairs (e.g. S1-S2, W1-W2, E1-E2). The visibilities were either 
lower than the visibilities in the other blocks on the same baseline or possessed an opposite 
or flat slope in comparison to other blocks on the same baseline. Data was judged to be dis¬ 
crepant via visual inspection and removed prior to either modeling or image reconstruction. 
At most, the rejected data only amounted to 1% of the total data for any epoch. 





2.2.3. 31 And 


The star 37 And was used as a calibrator, in addition to other calibrators during the 
same night, for the 2008, 2009 and 2010 observing runs. Che et ah (2012) discovered that 37 
And is a high contrast binary based on sinusoidal-like variations in the closure phase with a 
few degree amplitude. 


The orbit of the 37 And companion has been fully characterized by Roettenbacher et ah 
(2015, in preparation). To study the effect this companion has on A And data, the visibility 
and closure phase of a single star (the assumed calibrator) and a binary system (the actual 
calibrator) are compared to each other. The single star is taken to be a uniform disk with 
an angular diameter of 0.789 mas. The binary system is composed of the single star and 
a point source with a flux ratio of 80 (primary/secondary) and a separation of 46.66 mas, 
the maximum component separation. It is very likely the separation was shorter during the 
observations in question, however the full separation is used to determine the large possible 
error introduced. The [u,v] coverage for this analysis is the same as Sep 2"''^+3^'^, which 
represents the densest coverage of any single epoch. 

Fig shows both the visibilities expected from the uniform disk and the binary. Also 
shown is the uniform disk visibility divided by the binary visibility as a function of base¬ 
line. The error introduced increases with baseline to a maximum of 4% with an average of 
~1.5% Given the typical multiplicative errors in squared visibilities are 10-15%, this effect 
is considered negligible. It should also be noted that this represents the worst contribution 
made by the binary and does not address the contribution of additional proper calibrators 
that will mitigate the above effect. 

Fig shows the expected zero closure phase for the single uniform disk. The red dia¬ 
monds indicate the closure phase signal produced by the binary. The standard deviation of 
this signal is 1.27°. Given the noise floor for A And measured closure phase is 1°, the effect 
of the binary is negligible. 


2.3. Photometric Observations 

Photometric observations of A And were obtained using the 0.4-m automated telescope 
at Fairborn Observatory operated by Tennessee State University. A And was observed 376 
times for 3.4 years ranging from Sep 20*^, 2008 to Feb 1'^*, 2012. Fig shows the measured 
time series. The differential time series was corrected for atmospheric extinction and trans¬ 
formed into the Johnson U hlter system. T And (U: 4.982) and k And (U: 4.137) were 
used as the comparison and check star, respectively. The typical photometric error is 6.0 
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millimag. These errors are computed from the standard deviations in the time series of the 
check star. 


2.3.1. A And Optical Light Curve 


The qualitative behavior of the time series (i.e. overall trends, changing variability 
amplitudes) is consistent with A And’s 11.1 yr stellar cycle (Henry et al. 1995). The time 
series can be split into 4 “seasons”; each season corresponds to photometric observations 
taken just after the interferometric observations of the same year. The offset between the 
beginning of the interferometric and photometric observations is ~1 A And rotation. The 
seasons are labeled for the year in which the observations were taken. The time series period 
and variability amplitude in each season is measured for comparison with the modeled light 
curves generated from the interferometric images. Table contains the date ranges of each 
season, along with the identified period and the peak-to-trough A V amplitude. 

The period of photometric variability for the time series in each season is identihed 
using the Plavchan-Parks algorithm (Parks et al. 2014[ ). Fig displays the light curve for 
each season folded to the most signihcant period. Uncertainties in the period are set by 
the width of a Gaussian £t to the most signihcant peak in the power spectrum. The four 
determined periods are 54.27 ± 0.032, 55.15 ± 0.91, 53.3 ± 1.1, and 53.3 ± 1.9 days. The 
doubled appearance in season 4 could arise from starspots primarily on longitudes separated 
by 180°. In the previous seasons, one or more starspots only dominate a single hemisphere. 
The average rotation period of A And is 54.02 ± 0.88 days, where the reported error is the 
standard deviation of the four measured periods. 


3. STARSPOT ANALYSIS 

Two different techniques are employed to characterize starspots from the observed in¬ 
terferometric data: a cool spotted stellar surface model and image reconstruction. These 
methods are independent of each other; the results of one technique were not used as a 
starting condition for the other. 


3.1. Spotted Star Model 

The stellar surface model is computed using an IDL code capable of modeling any 
number of circular cool or hot starspots on a limb-darkened surface. In addition, the code 
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accounts for the effects of fore-shortening on starspots located away from the sub-stellar 
point. The free parameters are the stellar angular diameter [6), power law limb-darkening 
coefficient (a), starspot covering factor (0), starspot latitude {b), starspot longitude (/), and 
starspot intensity ratio (/). The code extracts model interferometric data by computing the 
Fourier transform of an artihcially generated stellar surface projected on the plane of the sky. 
The sampling for the Fourier transform is taken from the [u,v] coverage of the observed data 
being modeled. The goodness-of-£t parameter is the equally weighted reduced average 
between the observed and modeled visibilities, closure phases and triple amplitudes. 

Changes in the angular diameter and, to a lesser extent, the limb-darkening coefficient 
have a large effect on modeled visibilities at spatial scales smaller than the hrst zero of the 
visibility function. As the starspot information is contained on these small spatial scales, 
accurately determining the stellar properties prior to searching for the starspot properties is 
needed. This is done by first combining all the interferometric data from 2010 and 2011 into 
a single OIFITS hie. The 6 and a are measured by modeling only the visibility data using a 
grid search with a parameter step size of 1 x 10“^. The parameter uncertainty is determined 
by the lowest reduced -|- 1 ellipse in the reduced space. The furthest points of the 
ellipse along each axis are taken as the error in 6 and a, respectively. Closure phases are not 
considered at this stage; A And is likely not rotationally (vsinz = 6.5 km/s) or Roche lobe 
distorted. 


Once 6 and a are known, model solutions for each epoch are computed using a multi¬ 


parameter downhill simplex minimization (Press et ah 1992, AMOEBA). Model solutions 
with one, two, and three starspots are run with the preferred model yielding the lowest 
y^ statistic. A fourth starspot model is only investigated if the presence of the additional 
starspot can be inferred from with the prior and subsequent epochs based on rotational 
grounds. Only cool starspots are modeled as these are the type to persist on time-scales 
of a stellar rotation. Early attempts with the AMOEBA algorithm on starspot models 
demonstrated that the solutions are biased by initial parameters and search scales. This is 
indicative of numerous local minima in the reduced y^ space along with the deeper global 
minimum. The search scale employed is roughly 10% of the physical range for each param¬ 
eter. For example, the range in allowable intensity ratios is from 0.5 to 1.0 so the search 
scale is set to 0.05. AMOEBA is only very effective at hnding an accurate solution once the 
search occurs in the global minimum. Therefore, a genetic algorithm is employed prior to 
running AMOEBA in order to hnd the approximate location of the reduced y^ minimum. 

A genetic algorithm (GA) is an iterative process through which a best solution is found 


by “evolving” an initial set of randomly chosen model solutions (members) (Charbonneau 


1995). The htness, or chance it will be used in the subsequent iteration, of each member is 
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determined based on the member’s reduced The “survival” process is determined via a 
roulette wheel scheme. The wheel is spun a number of times equal to the population size. 
The probability the wheel will choose a member to survive is proportional to the member’s 
htness. Therefore the next population will, in theory, be composed of model solutions with 
lower on average. This new population is “evolved” via two different random methods: 
crossover and mutation. Crossover takes sections of a parameter value and swaps it with 
another parameter value. For example, solution A has a latitude of 45.12° and solution 
B has a longitude of 12.57°. Crossover can swap the digits after the decimal to yield a 
new latitude of 45.57° and longitude of 12.12°. Mutation causes a section of the parameter 
value to change randomly. Using the previous example, the latitude 42.12° could mutate 
to become 49.12°. Both crossover and mutation are applied with a frequency of 90% and 
1%, respectively. The htness of the new population is determined and the entire process is 
iterated until the average htness drops below a convergence criterion. 

The reason AMOEBA is still needed after the GA is the GA’s inability to converge to 
the exact minima in x^ space. During the GA and AMOEBA implementations, both the 
stellar angular diameter and limb-darkening coefficient are kept hxed. 

Errors to model starspot parameters are found by randomly varying each parameter of 
the best-ht model and then re-running the AMOEBA search algorithm. All the parameters 
are uniformly randomized simultaneously within a given range. 0 and / are allowed to vary 
by at most ±0.1 while b and I are varied by 1.8°. This procedure, which is meant to explore 
the reduced X 2 -space around the best-ht solution, is performed ten times. The parameter 
errors correspond to the standard deviation of the ten values yielded by the above procedure. 
If a trial solution has a x^ better than the initial hnal solution, this trial solution becomes 
the hnal solution. 

Since the observations occur in the Rayleigh-Jeans tail of the spectral energy distri¬ 
bution, the temperature ratio between a starspot and the photosphere should be a linear 
function of the corresponding intensity ratio. Therefore, the intensity ratio and its error are 
directly translated into the starspot temperature ratio, T/j, and corresponding error. 


3.2. SQUEEZE: Image Reconstruction 

Model-based imaging is a very ehective tool in determining starspot properties; however 
it is limited by the assumptions used to create the model (e.g. circular starspots). Model- 
independent imaging, or image reconstruction, on the other hand has freedom to portray 
more realistic starspot shapes and sizes. The main hurdle faced by image reconstruction is 
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an ill-posed inverse problem bronght on by working in a data starved regime, incomplete 
[n,v] coverage, and atmosphere-corrupted Fourier phase information. Image reconstruction 
algorithms are all based on a regularized maximum likelihood paradigm that reconciles a 
statistic with k number of regularization statistics, or regularizer (i?k), modulated by a 
user-dehned weighting parameter (/ik)- 


i = argmax 
ieM" 


K 


X 


(i) -|- ^ ^ k'k^k (^) 


k=l 


( 1 ) 


Image reconstructions abide by two restrictions: the flux in a particular pixel must be non¬ 
negative and the flux of the reconstructed image is normalized to unity. 


The image reconstruction code SQUEEZE]^ is used on the A And data sets (Baron et ah 


2010); this code is an evolved version of the reconstruction program MACIM (Ireland et ah 


2006). The SQUEEZE initial state is a 50 x 50 pixel array of a 2.777 mas (the measured 6 for 


A And) uniform disk containing 4000 flux elements. After the initial state is set, SQUEEZE 
uses a multi-thread approach with each thread Ending the best image by moving flux from 
pixel to pixel iteratively via a simulated annealing algorithm. The total variation regularizer 


is chosen as it is designed to minimize brightness gradients across the surface (Rudin et ah 


1992). This favors a conservative surface image with a few large starspots as opposed to 


many smaller ones. Experiments with other regularizers such as Laplacian regularization 
and the io sparsity norm were attempted but failed to produce images with reduced lower 
than total variation alone. This is consistent with the findings of Renard et ah (2011) who 
found total variation outperformed four other regularization methods. 


A final image reconstruction is the average of ten images generated by SQUEEZE. 
The initial conditions in regards to the initial image (uniform disk), choice of regularizer 
(total variation), and hyper-parameter are identical for each ten images. The random seed 
used for the Markov chains to create the image is unique for each image. Prior to being 
averaged, the ten images were registered using SPLASHj^to account for any changes in the 
image photocenter. The averaging of ten images attempts to minimize the effect of artifacts 
potentially caused by the reconstruction process. Starspot parameters are extracted by 
fitting a circular aperture over identified starspots. This process is illustrated in Fig[^ The 
aperture size provides the covering factor and the location of the aperture center provides 
the starspot latitude and longitude. As the starspot edge is difficult to quantify and the 


^https: //github.com/fabienbaron/squeeze 
^https: //github.com/fabienbaron/splash 
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starspot may be irregular in shape, 0 is considered a lower bound. The intensity ratio is 
calculated by dividing the flux at the aperture center with a flux measurement of the “quiet” 
photosphere. The quiet photosphere is identified as a part of the stellar surface devoid of 
flux gradients. The circular aperture is fit to the reconstructed starspots by eye. 

In addition to the creation of a final averaged reconstructed image, an image representing 
the standard deviation of the ten iterations is created. The detection strength, ads, of the 
starspot is computed using this mean standard deviation. A circular aperture is placed on 
the quiet photosphere with a size equal to the minimum angular resolution. The detection 
strength is the mean flux within this aperture subtracted from the mean flux within the 
starspot aperture and then divided by the standard deviation. 


3.2.1. Image Artifacts 


A large problem facing image reconstruction is the believability of surface features. 
Certain results from imaging can be easily dismissed as artifacts such as hexagonal stellar 
surfaces or a repeating symmetric patterns of bright spots that is inconsistent with a rotating 
surface. Certain results can be accepted as believable such as dark starspots common between 
the reconstructed and modeled images due to the independence between the two imaging 
methods. However, since single bright starspots were not modeled, this test is unavailable. 
Therefore a method must be devised to determine if any imaged bright starspots are real (e.g. 
flares, plages) or a consequence of the data sampling and/or image reconstruction process. 


The method used in this analysis is to create synthetic reconstructed images based on 
the model images. This is done by extracting visibilities, triple amplitudes, and closure 
phases from each model image with identical data sampling and SNR as the measured data 
using the OIFITS-SIM tool^ The same reconstruction procedure described in § 3^ is used 
to create synthetic reconstructed images from these data. Artifacts due to miscalibrated 
observables will be features seen in the reconstructed image, but are absent in both the 
simulated and the model images. Artifacts due to the [u,v] coverage and the reconstruction 
process will be features seen in both the reconstructed and simulated images, but not in the 
true image. 


^https://github.com/bkloppenborg/oifits-sim 
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DISCUSSION 


4.1. A Andromedae Stellar Properties 


The angular diameter and limb-darkening coefficient are determined via the modeling 
The initial value of 6 for the AMOEBA code is set to 2.75 mas as 


described in 5 3.1 


determined from the A And ( V-Kg) color (van Belle|1999 ). An initial a is found by matching a 


power-law £t to a four parameter £t from Claret & Bloemen (2011) given the coefficients for a 


star with Tes = 4750 K and log(g) = 3.0 dex (cgs). This yielded a result of a = 0.22 consistent 


with results from other power law hts to interferometric data of late type giants (Wittkowski 


et al.||2002 2006). The search scales were roughly 10% of the initial values. The hnal results 
are 6 = 2.759 ± 0.050 mas and a = 0.229 ± 0.111. The uncertainty is measured from the 
lowest + 1 ellipse in the reduced space. The error bars correspond to the furtherest 
points of the ellipse along each axis, respectively. Figj^shows the reduced space for 6 and 
a with the white cross marking the best solution. A slight degeneracy exists between these 
two parameters as seen by the elliptically-shaped gradients. At a Hipparcos trigonometric 


distance of 37.87 ± 0.21 mas (van Leeuwen 2007), this angular diameter corresponds to a 


linear radius of 7.831]'^g;Qg5 R, 


'©• 


4.2. A Andromedae Starspot Properties: 2010 Data Set 


Between Aug 2”'^ and Sep 11*^, 2010, 6 epochs of data were obtained via the strategy 
The number of [u,v] points obtained for each of the two night combined 


described in 5 2.2 


5 epochs ranged from 624 to 1128 with the densest converge obtained by the combination 
of Sep 2’^'^ and 3^*^ (see Table [^. The number of [u,v] points for the single night epoch 
on Sep 10*^ is 336. Fig 0 and Fig show the distribution of [u,v] coverage obtained for 
each pair of observations. The six epochs are spaced with a cadence between six to nine 
days corresponding to 10.9% to 16.4% of the measured rotation period; signihcant apparent 
starspot motion is expected between epochs. The complete observing run spans 71% of one 
complete A And rotation period. 


Fig shows a distinct non-zero closure phase signature across most sampled spatial 
scales pointing to the existence of surface asymmetries. This signature is present in all 
six epochs. The measured closure phase distribution differs between epochs supporting the 
hypothesis that the apparent starspot conhguration evolves as the star rotates. In addition, 
an unspotted model image does not £t the interferometric data for each epoch in 2010 with 
the reduced ranging between 15.4 and 44.0. 















15 


The best-fit models for each epoch contain between two to fonr cool starspots. Fig 10 


contain the model, reconstructed and simulated images for each epoch. The model reduced 
range between 2.97 to 6.02 with the best £t occurring for Aug 10*^ and 11*^. As an 
ensemble, the covering factor for individual starspots, 0, ranges from 4.0 to 21.8% with a 
median value of 7.6%. The covering factor errors range 2.6 to 5.5% with a median error of 
4.1%. The temperature ratio, ranges from 0.756 to 0.926 with a median value of 0.853. 
The temperature ratio errors range from 0.007 to 0.042 with a median error of 0.015. The 
errors in both latitude and longitude are nearly identical and range from 0.75 to 7.8° with a 
median error of 1.9°. 


The reduced for oach reconstructed image is < 1.01. Good qualitative agreement 
exists between the model and reconstructed images (see Fig 10). Table contains the 
measured starspot parameters from both the model and reconstructed images. This table 
also contains the difference between the parameters for the same starspot in the model and 
reconstructed images. 


The Epoch 6 model image contains two cool starspots. Although both starspots are 
listed in Table the starspot located near the northeastern limb {h\ 17.5°, l\ -40.9°) is 
excluded when discussing ensemble starspot properties. The rationale for this decision is: 
the starspot is nearly twice the size (0 = 44%) of the next largest modeled starspot, it is 
the warmest starspot (Tj:{ = 0.925), and is not conhdently seen in the reconstructed image. 
This “starspot” could be a widely spread patch of starspots with the covering factor of each 
individual starspot falling below the resolution limit . The model attempts to reconcile the 
interferometric signature these starspots produce through the addition of this larger, warmer 
starspot. 


The reconstruction process does not identify starspots on the stellar limb nearly as well 
as the model. For epochs with visual agreement between model and reconstructed starspots, 
0 and T/j estimates are within one error bar in nearly every case. The agreement suffers for 
position estimates with starspots in reconstructions consistently further east and/or south 
than their model image counterparts. 


The simulated images and the observing cadence are used to help identify potential 
artifacts in the reconstructed images. In Epochs 2 through 6, a number of warm starspots 
are observed in the reconstructed image evenly spaced around the stellar limb. These are 
rejected as artifacts due to their symmetry and constant position despite the rotating surface. 
The origin of these artifacts may be due to the [u,v] sampling since the pattern of the warm 
starspots is similar to the pattern of tightly clustered points in the [u,v] plane (see Fig|^. The 
southern cool starspot in Epoch 2 is rejected as the covering factor is below the resolution 
limit. The central warm starspot in the reconstructed images of Epochs 3 and 4 are rejected 
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as artifacts; they are reproduced in the simulated images and are contrary to a rotating 
surface. The brighter southern pole in Epoch 6 is reproduced in the simulated image and 
thus rejected. The warm southwestern starspot Epoch 1, however, cannot be rejected as 
false as it is not reproduced. 


4-2.1. Comparison with the 2010 light curve 


Determining if the starspot imaging provides a consistent picture with the measured 
photometric time series will provide further conhdence in the interferometric results. Season 
2010 spans 121.8 days and has peak-to-trough photometric AV amplitude of 0.099 mag. The 
time series begins 10 days after the interferometric observations providing no simultaneous 
coverage. The seasons spans ~2.2 rotation periods, and varies from one rotation to the 
next indicating starspot evolution between rotations. Also the light curve seems bimodal 
suggesting starspots on two active longitudes (0 = 0.0 and 0.4). Active longitudes, or 
longitudes of preferential starspot formation, have been associated with magnetically active 
stars ( Berdyugina|2005 ). This could be an example of active longitude migration with respect 


to the star’s rotational frame of reference (Jetsu et ah 1993 Lanza et ah 1998 Berdyugina 


& Tuominen 1998) 


A modeled light curve can be constructed from the best-£t models for each epoch. 
The change in intensity between an unspotted star and the modeled surface is measured, 
converted into a A magnitude and then scaled for comparison with the observed time series. 
The scaling is done through an additive constant that shifts the modeled time series to 
the approximate values of the photometric time series. A multiplicative constant is used 
to expand the variability amplitude to be comparable to the photometric time series. A 


modeled light curve is computed and shown in Fig overlaid on top of the photometric 
time series folded to the 53.35 day period. There is good visible agreement in all epochs 
except of Epoch 5, which is brighter than expected. It is difficult to explain the discrepancy 
based on [u,v] coverage since this epoch had the densest coverage. Data quality does not 
seem to be a viable explanation as the errors are not signihcantly larger than other epochs 
and the model reduced is one of the lowest all six epochs. 


4-2.2. Tracing Rotation in the 2010 Data Set 

Multi-epoch starspot imaging has the potential to trace the stellar rotation via apparent 
starspot motion. If this motion can be observed, a rotation period can be measured and 
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compared to the period derived from photometry. In addition, the stellar rotation axis can 
be fully described in both inclination and position angle. Neither Doppler imaging or light 
curve inversion has the capacity to determine the latter two quantities and, in fact, the 
inclination angle must be assumed in both cases. 


The observing baseline for the 2010 data set spans 71% of the photometrically deter¬ 
mined rotation period. The cadence will carry starspots ~47° across the stellar surface 
between epochs assuming a negligible amount of differential rotation. 


Individual starspots identihed in different epochs are labeled A through (7 in Fig [T^ 
Identihcation of the same starspot in different epochs is judged via visual inspection. Four 
starspots (C, D, E, F) are believed to be seen in three epochs and therefore provide the 
most useful constraints on the rotation and inclination angle. Changes in covering factor 
and temperature ratio are described in Table along with the computed rotation period 
based on the measured angle between starspot positions from one epoch to the next. The 
change in 0 and should be small since these starspot characteristics are not expected to 
evolve during one stellar rotation. The median A0 and A Tji are 0.5% and 0.014, respectively. 
As these are below the median (j) error of 4.6% and median Tji of 0.016, it is safe to claim 
that these starspots are not signihcantly evolving. The range in measured rotation periods is 
44.3 to 78.7 days with an average rotation period of 60 ± 13 days; the error is the standard 
deviation of the individual measured rotation periods. The starspot determined period is 
consistent with the photometric rotation period of 54.02 ± 0.88 days due to the large error 
bar. The large spread in starspot determined rotation periods is likely a consequence of four 
starspots observed only three times over an incomplete rotation. 


Fig shows each of the starspots plotted by longitude vs. latitude overlaid by ellipse 
hts. Estimations of the inclination and position angle of the rotation angle can be made 
by measuring these elliptical paths. A starspot being carried across the stellar surface via 
rotation will appear to travel along an elliptical path when viewed in two dimensions. The 
position angle, T, is simply the tilt of this ellipse counterclockwise from north (up). The 
inclination angle, i, is the inverse sine of the ellipse eccentricity. If the star is viewed face-on, 
then the starspot will appear to traverse a circular path. Conversely if the star is viewed 
edge-on, the starspot will appear to traverse a line. An ellipse is £t via visual inspection to 
starspots C through F as these spots have been identihed in at least 3 epochs. The average 
T and i are 19 ± 8.1° and 75 ± 5.0°, respectively. The errors are the standard deviations 
of measured values. This inclination angle is higher than the 60° assumed by F08, but is 
consistent within the uncertainties in the previous inclination estimate 601^5° by D95. The 
rotation axis is measured to be coming out of the plane of the sky in the northern hemisphere. 
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4.3. A Andromedae Starspot Properties: 2011 Data Set 

Between Sep 2’*'^ and Sep 24*^, 2011, six epochs of data were obtained for A And nsing 
all six telescopes simnltaneonsly. In order to maximize the [n,v] coverage, observations were 
taken for as long as delay lines were available. Only one night of data was acqnired per epoch 
as opposed to the combined days in 2010. Table [^contains the nnmber of [n,v] obtains points 
per observation. The nnmber of [n,v] points achieved ranges from 200 to 864 with the densest 
coverage obtained on Sep 14*^ (see Fig 14). The observing cadence was fonr or hve days 


corresponding to 7.3% and 9.2% of the rotation period; all six observations span 40.4% of 
one rotation period. 


Fig shows a distinct non-zero closnre phase signatnre across most sampled spatial 
scales in all six epochs pointing to an evolving conhgnration of snrface asymmetries. In 
addition, an nnspotted model image does not £t the interferometric data for each epoch in 
2011 with the rednced ranging between 9.3 to 16.7. 

The best-£t models for each epoch contain one or two cool starspots. The model rednced 
range between 3.14 to 6.57 for these epochs with the best £t occnrring for Sep 2"^^. The 
starspot properties are listed in Table Fig IT contain the model, reconstrncted and 
simnlated images for each epoch. As an ensemble, the covering factor, 0, ranges from 10% 
to 17% with a median valne of 12%. The errors in these valnes range from 2.5% to 6.4% 
with a median error of 5.3%. The temperatnre ratio, Tr, ranges from 0.797 to 0.867 with 
a median valne of 0.842. The temperatnre ratio errors range from 0.007 to 0.016 with a 
median error of 0.014. The errors in both latitnde and longitnde are nearly identical and 
range from 0.63 to 5.7° with a median error of 1.4°. 

The rednced for each of the reconstrncted images is at or below 1.01. Good qnalitative 
agreement exists between the model and reconstrncted images, however some exceptions do 


exist (see Fig 16). It is nnclear why the western starspot seen in model images is not recovered 
in the reconstrncted images for Sep 6*^ and 10*^. Nor is it clear why there is poor agreement 
between the Sep 19*^ model and reconstrncted images. In general, the agreement between 
starspot properties is not as good as it is in the 2010 data set. The reconstructed covering 
factor is always smaller than the modeled covering factor in each epoch, however, this is not 
surprising as this parameter should be considered a lower bound. In an opposite trend than 
2010, the reconstructed b is further north where it deviates from the model image. No trend 
exists for when the agreement in I exceeds the uncertainties. 


The reconstructed images in Epochs 1, 3, and 5 contain a ring of warm starspots around 
the stellar limb. These starspots are rejected as artifacts caused by [u,v] sampling due to 
their symmetry and constant location between epochs, which is contrary to starspots on a 
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rotating surface. The following starspots are rejected as artifacts by the presence of similar 
features in the simulated image. In Epochs 2 and 4, the warm starspots in the northeast 
and southwest are rejected as artifacts. The warm northeast starspot in Epoch 3, the warm 
central starspot in Epoch 4, and the cool southern starspot in Epoch 5 are all rejected 
as artifacts. The non-circular stellar disk in Epoch 6 is most certainly an artifact due to 
the limited number of [u,v] points (200). In addition, the shape of this disk resembles the 
conhguration of [u,v] points (see Fig 14). 


4-3.1. Comparison with the 2011 light curve 

Again the measured photometric time series is compared to the imaging results to 
provide further conhdence in the interferometric results. Season 2011 spans 134.9 days and 
has a peak-to-trough AV magnitude of 0.057 mag. This season spans ~2.5 rotation periods 
and overlaps the hrst two epochs of interferometric observations. The folded light curve is 
double-peaked; this is potentially due to two active longitudes separated by ~180°. 

A modeled light curve is computed and shown in Fig [T^ overlaid on top of the photo¬ 
metric time series folded to the 53.30 day period. The modeled photometry in each epoch 
is given by the red asterisks and is agreement with the observed photometry in Epochs 2, 
3, and 5. The agreement for Epoch 1 is not as convincing, but not completely inconsistent. 
However, the modeled photometry for Epochs 4 and 6 visibly disagrees with the observed 
photometry. The hlled blue circles represent the modeled photometry, if the parameters for 
a single starspot are changed. The altered starspots are B in Epochs 1 and 4 and C in 
Epoch 6. The starspot covering factor in both Epochs 1 and 4 is lowered to 11.8%, which 
is identical to that of Epoch 3 where the starspot is most centrally located. Additionally 
the Epoch 4 starspot T/j is lowered to that of starspot B in Epoch 3 (0.797). The starspot 
covering factor in Epoch 6 is lowered to 10.4%, which is the same as the 0 for the same 
starspot in Epoch 5. 


4-3.2. Tracing Rotation in the 2011 Data Set 


The observing baseline for the 2011 data set spans ~41% of the photometrically de¬ 
termined rotation period. The cadence will carry starspots ~7% across the stellar surface 
between epochs assuming a negligible amount of differential rotation. Fig 18 shows a com¬ 
pelling pattern of stellar rotation by three starspots labeled A through C. Starspot B is 
seen in all 6 epochs and provides the best estimates of both the stellar rotation and rotation 
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axis. The observing strategy behind the 2011 data set was designed to provide an increased 
number of measurements for any individual transiting starspot (s). This is done to shrink 
the uncertainties in the estimates of the rotation axis computed from the 2010 data set. The 
uncertainties arose from a sparse number of measures per transiting starspot. 

Table contains the changes in 0 and T/j for starspots A, B, and C and the computed 
rotation period based on the measured angle between starspot positions from one epoch to 
the next. The largest change in cj) and for any starspot over the observed rotation is 
4.6% and 0.063, respectively. As these are comparable to the median cj) error of 5.3% and 
median Tn of 0.014, the starspots do not appear to be signihcantly evolving. The range in 
the measured rotation period based on apparent starspot motion is 43.4 to 64.6 days with 
an average rotation period of 54.0 ± 7.6 days, which is nearly identical to the photometric 
rotation period of 54.02 ± 0.88 days. The error is the standard deviation of the individual 
rotation periods. 


Fig shows each of the starspots plotted by latitude vs. longitude overlaid by ellipse 
fits. For starspots A and B, an ellipse is fit via visual inspection. The C starspot is excluded 
from this analysis due to having only two data points. The average T and i are 28 ± 1.2° and 
69 ± 1.4°, respectively. The rotation axis is tilted out of the plane of the sky in the northern 
hemisphere. These values and orientation are nearly identical to the estimates found by the 
2010 data sets. 


4.4. Comparing Results With Previous Work 

Having demonstrated that starspot properties can be measured for A And using inter¬ 
ferometric observations, these results are compared to the results of previous investigations. 
D95 created a surface map of A And via a matrix LCI technique using Johnson BV light 
curves spanning one rotation period. D95 models the observed light curve using 2 starspots 
with a Tr = 0.83 ± 0.06. One starspot is located at 6 = 50° with a 0 = 8%. The other 
starspot is located at 6 = 20° with a 0 = 4%. The starspots are separated by 140° in lon¬ 
gitude. Both the latitudes and covering factors for these starspots are consistent with those 
identified in this work. However, the temperature ratio in D95 is significantly less than that 
measured for both the 2010 (median Tr = 0.961) and 2011 (median Tr = 0.958) data sets. 

A more recent study of the starspot properties of A And was performed by F08. They 
use a 2-component LCI method using Johnson F-Band photometry coupled with spectral 
line depth ratios to create a map of starspots on A And. The results of F08 are very consistent 
with D95 with the modeled surface containing 2 cool starspots each with a Tr = 0 . 815 °g®Q 3 g. 
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The covering factors for the two starspots are 8.7% and 3.6% located at latitudes 57° and 9°, 
respectively. The starspots are separated by 81° in longitude. As with D95, the positions and 
covering factors of the starspots identified by F08 are consistent with the starspots imaged 
in this work. However, the temperature ratio is significantly less. One additional difference 
between F08 and this work, as well as D95, is the modeling of 2 plage regions by F08. These 
bright regions are similar in size to the modeled cool starspots. The plages are also in similar 
locations only offset to the starspots by ~20° in longitude and ~7° in latitude. By design 
no hot starspot/plage regions were modeled in this work, however the reconstructed images 
show no convincing evidence these exist in any observed epoch. 


5. SUMMARY 


A Andromedae, a bright (V: 3.872 mag, H: 1.501 mag) G8 giant, has a long recorded 
history of consistent sinusoidal-like photometric variability. This variability is believed to 
result from the rotational modulation of cool starspots. Using light curve inversion tech¬ 


niques, the presence of starspots have been indirectly revealed (Donati et al. 1995 Frasca 


et al.||2008 ). Long baseline optical/near infra-red interferometry has the potential to directly 


detect and measure the characteristics of these starspots. LBI has already directly imaged 
a number of astrophysical systems (e.g. close binaries, circumstellar disks, rapidly rotating 
stellar surfaces) with unprecedented angular resolution. 


A And was observed using the MIRC beam combiner on the CHARA Array for 26 days 
spanning from Aug 17*^, 2008 to Sep 24*^, 2011. The observing strategy evolved over time 
due to upgrades in the MIRC beam combiner. Contemporaneous photometric observations 
are also available from Sep 20*^, 2008 to Jan 20*^, 2011. The photometry provides an 
independent relative estimate of starspot coverage that can be compared to the imaging 
results. 


The CHARA array observations in 2010 employed an observing strategy designed to 
maximize the [u,v] coverage. The images created from the 2010 data feature one to four 
starspots at each of the 6 epochs. As an ensemble, the median value in the starspot cov¬ 
ering factor is 7.6%. The median value of the temperature ratio between starspot and the 
surrounding photosphere is 0.853. The starspot properties extracted from the reconstructed 
images are consistent with the modeled results to within the error bars. The observing ca¬ 
dence between the 6 epochs in 2010 is between 6 to 9 days corresponding to 10.9% to 16.4% 
of the rotation period. The observations span 71% of one rotation cycle. By measuring 
the apparent motion of four starspots for 3 epochs, the rotation period is measured and the 
rotation axis is determined. The rotation period based on starspot motion is 60 ± 13 days. 
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which is consistent with the photometric rotation period of 54.8 ± 1.9 days. The rotation 
axis is tilted ont of the plane of the sky in the north with an inclination of 78.0 ± 1.5° and 
a position angle of 20.0 ± 6.8° 

A photometric F-Band time series is available beginning 10 days after the last interfero¬ 
metric observation in 2010. The short-cadence time series spans approximately two rotations 
of A And. The flux variability in the modeled time series follows the behavior of the photo¬ 
metric time series quite well when the proper scaling factors are applied with the exception 
of a single epoch. The discrepancy in this epoch is difficult to explain based on [u,v] coverage 
or data quality. 

In 2011, A And was observed for a single night on 6 different nights using 6 telescopes 
simultaneously. This strategy substantially increased the number of visibilities and closure 
phases obtained for each block of data compared to 2010. However, the [u,v] coverage 
decreased by appropriately a factor of 2 since the data were only composed of a single night 
per epoch. Consequently, the model and reconstructed images are not as consistent as those 
created from the 2010 data. One to two starspots are identified in the model images for each 
epoch. As an ensemble, the median value of the starspot coverage is 12% and the median 
temperature ratio is 0.842. The 2011 observing cadence is 4 or 5 days corresponding to 7.3% 
and 9.2% of the rotation period spanning ~40% of one rotation. One starspot is imaged in 
all of the 6 epochs and is used to measure a rotation period of 54.0 ±7.6 days, which is 
nearly identical to the photometric rotation period of 54.8 ± 1.9 days. The rotation axis 
points out of the plane of the sky in the north with an inclination of 77.98 ± 0.18° and a 
position angle of 23.0 ± 6.4°. These values are in agreement with the orientation computed 
from the 2010 data set. 


The photometric F-Band time series starting in 2011 overlaps the last two interfero¬ 
metric epochs and spans ~2.5 rotations of A And. The flux variability in the modeled time 
series follows the behavior of the photometric time series only for half of the epochs. The 
discrepant epochs can be made to agree with the photometry if the covering factor and/or 
the temperature ratios are changed for a single starspot. The parameters are changed to 
be identical to parameters identified for the same starspot in a different epoch where the 
starspot is more centrally located. 


Using the combined 2010 and 2011 data sets the measured angular diameter of A And 
is 2.759 ± 0.050 mas with a power law limb darkening coefficient of 0.229 ± 0.111. At a 
Hipparcos trigonometric distance of 37.87 ± 0.21 mas, this yields a linear radius of 7.83l/^o;og5 
R©. The angular diameter is consistent with that previously found by Nordgren et ah (1999) 
(2.66 ± 0.08 mas). 





The 2010 and 2011 data sets have provided a convincing picture of starspots on the sur¬ 
face of A And in support of the closure phase information and the variable light curve. The 
agreement between the modeled and reconstructed starspot properties is typically within 
one error bar. In addition, the starspots produce a flux variability which is consistent with 
that observed photometrically just subsequent to the interferometric observations. There 
is evidence to suggest that starspots imaged in one epoch are again imaged in subsequent 
epochs. This provides starspot determined rotation periods consistent with the photometri¬ 
cally identihed periods. This method, additionally, allows for the complete characterization 
of the rotation axis of A And. 
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Fig. 1.— Top-. Shown here are both the visibilities expected from a uniform disk (0.798 
mas) and the binary (point source secondary, flux contrast = 80, separation = 46.66 mas). 
bottom,-. This contains the uniform disk visibility divided by the binary visibility as a function 
of baseline. 
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Fig. 2.— Shown here are both the closure phases expected from a uniform disk and the 
binary. As expected the closure phases for the uniform disk are zero, while the standard 
deviation of the binary closure phases is 1.27°. 
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Fig. 3.— Time series photometry for A And ranging from Sep 2008 to Feb 2012. The tick 
marks along the bottom axis represent the interferometric observations. The brightening 
trend of the time-series is due to the 11.1 yr stellar cycle (Henry et ah 1995) 
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Fig. 4.— Top Left: Season 2008 time series folded to a period of 54.27 ± 0.032 days. Top 
Right: Season 2009 time series folded to a period of 55.15 ± 0.91 days. Bottom Left: Season 
2010 time series folded to a period of 53.4 ± 1.1 days. Bottom Right: Season 2011 time 
series folded to a period of 53.3 ± 1.9 days. 
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Fig. 5.— Shown is a closeup of the SQUEEZE reconstruction for the Sep 2*^, 2011 data 
near an apparent starspot. The black circle on the right shows the aperture used to ex¬ 
tract starspot properties from reconstructed image. The black circle on the left shows the 
aperture over the “quiet” photosphere. The “quiet” photosphere is dehned as a part of the 
stellar surface devoid of flux gradients. The size of the aperture is identical to the minimum 
achievable angular resolution. 
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Fig. 6.— The best fit 6 and a are determined throngh a grid search with a parameter step 
size of 1 X 10“^. The rednced space for 9 and a is shown with the white cross marking 
the position of best £t. A slight degeneracy exists between these two parameters as seen by 
the elliptically-shaped gradients. 
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Fig. 7.— The [u,v] coverage obtained for the first three epochs in 2010. A - Aug 2’^'^; B 
Aug 3’’'^; C - Aug 10*^; D - Aug 11*^; E - Aug 18*^; F - Aug 19*^ 
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ig. 8.— The [u,v] coverage obtained for the second three epochs in 2010. A - Aug 24*^; B 
Aug 25*^^; C - Sep 2’^'^; D - Aug 25*^^; E - Sep lO*'* 
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Fig. 9.— The observed closure phases for the 2010 data set. The errors bars have been 
excluded for clarity. The distinct non-zero closure phase signature points to surface asym¬ 
metries. The differences in the closure phase between nights indicates an evolving asymmetric 
surface pattern from night to night. 
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Fig. 10.— Stellar surface images for the 2010 data set. The top row contains the model 
images, the middle row contains the reconstructed images, and the bottom row contains 
the simulated images. The white dot in the lower right hand corner represents the 0.4 mas 
resolution limit for the CHARA array. 
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Fig. 11.— The gray diamonds correspond to the F-Band time series phased to a period of 
53.3 ± 1.9 days. The red asterisks represent the photometry taken from the best-£t models 
for the 6 epochs. The numbers along the bottom axis indicate the epoch. 
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Fig. 12.— The best fit models for each epoch in 2010. In each model, the starspot(s) are 
labeled {A through G) to indicate the same starspot as seen in each rotational phase. 
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Fig. 13.— Ellipse fits to starspot positions in the 2010 data sets. The dotted line corresponds 
to spot C. The solid line corresponds to spot D. The dot dashed line corresponds to spot 
E. The dashed line corresponds to spot F. The average computed position angle, <h, and 
inclination angle, f, from these hts are 19 ± 8.1° and 75 ± 5.0°, respectively. The solid line 
circle corresponds to the circumference of A And. 



V (*10® rad"') 








Closure Phase (deg) 


- 38 - 


200 


100 


0 


-100 


-200 



1 . 5 * 10 ® 2 . 0 * 10 ® 2 . 5 * 10 ® 3 . 0 * 10 ® 3 . 5 * 10 ® 4 . 0 * 10 ® 

Spatial Frequency (rad"") 


Fig. 15.— The observed closure phases for the 2011 data set. The error bars have been 


removed for clarity. The distinct non-zero closure phase signature points to surface asym¬ 
metries. The differences the closure phase between nights indicates an evolving asymmetric 
surface pattern from night to night. 
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Fig. 16.— Stellar surface images for the 2011 data set. The top row contains the model 
images, the middle row contains the reconstructed images, and the bottom row contains 
the simulated images. The white dot in the lower right hand corner represents the 0.4 mas 
resolution limit for the CHARA array. 
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Fig. 17.— The gray diamonds correspond to the F-Band time series phased to a period of 
53.3 ± 1.9 days. The red asterisks represent the photometry taken from the best-£t models 
for the 6 epochs. The hlled blue circles represent model photometry where the values of </> 
and/or T^j have been altered (see text). 
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Fig. 18.— The best fit models for each night in 2011. In each model, the starspot(s) are 
labeled (T, B, and C) to indicate the same starspot as seen in each rotational phase. 
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APPENDIX. THE EARLY IMAGING ATTEMPTS 

This appendix discusses the failed attempts at high fidelity surface imaging to illustrate 
the effect of insufficient [u,v] coverage on interferometric starspot imaging. Any successful 
program will need to be designed to maximize this coverage as much as possible. The 
inconsistent results measured in 2008 are first discussed and then followed by a discussion of 
the intriguing, yet inconclusive results measured in 2009. 


A-1. THE 2008 DATA SET 


In 2008, two observing runs of A And were performed with one in August and the 
other in September. Both runs employed “snapshot” observations using the S1-E1-W1-W2 
telescopes. These observations consist of only two or three bracketed observations per night. 
The August run was composed of observations taken on five consecutive nights between 
the 17*^ and the 21^*. The [u,v] coverage achieved ranged from 48 to 144 data points 
with the densest coverage obtained on Aug 18*^. Fig A-1 contains the plots of these [u,v] 
configurations. The September run was composed of two observations taken a week apart 


on the 20*^ and the 27*^. The [u,v] coverage achieved was 72 points for each night. Fig A-1 
contains the [u,v] configurations for these nights. 


Fig |A-2| shows a distinct non-zero closure phase signature across most sampled spatial 
scales for all five days. Unlike in 2010 and 2011, the closure phase signature changes only 
slightly from night to night, consistent with the hypothesis of surface asymmetries that have 
not evolved on short time scales. An unspotted model image yields an extremely poor fit to 
the interferometric data for each epoch with the reduced ranging between 5.6 to 18. 


The reduced for these models are all below 2.85 with the lowest reduced y^ (1-14) 
occurring on Aug 19*^. Despite this, these images do not present a consistent starspot 
configuration as can be seen in Fig |A-3[ This figure contains the model, reconstructed, and 
simulated images for Aug 17*^, 18*^, and 19*^. The one day cadence should present nearly 
identical surface images as the surface only rotates 1.8% from night to night. Starspot 
evolution on this time scale is not typical for magnetically active stars (Berdyugina| 2005 


Strassmeier[|2009). Additionally, the reconstructed images do not even contain any conclusive 


evidence for starspots. 


The likely factors contributing to the non-detection of a consistent starspot configuration 
include, but are not limited to, poor [u,v] sampling, mis calibration and the influence of the 
binary 37 And used as a calibrator. 
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Fig. 19.— Ellipse fits to starspot positions in the 2011 data sets. The dashed line corresponds 
to the £t to spot A. The dash dot dot line corresponds to spot B. The average computed 
position angle, <h, and inclination angle, i, from these fits are 28 ± 1.2° and 69 ± 1.4°, 
respectively. The solid line circle corresponds to the circumference of A And. 


44 


The phased photometric time series (Fig|^ indicates these interferometric observations 
where taken ~11 days after maximum brightness. As shown in the study by [Henry et ah 


(1995), the maximum brightness of A And can vary by as much as 0.05 mag in the F-Band. 
This would indicate that even at time of maximum brightness the visible surface still contains 
starspots. Therefore, the inconsistent images are probably not due to a insufficient starspot 
presence. 


Thus, despite strong evidence for starspots on the surface of A And during these epochs, 
from both measured non-zero closure phases and the variable light curve, the [u,v] coverage 
using 4 telescopes on a single night is insufficient to confidently determine starspot properties. 


A-2. THE 2009 DATA SET 


The A And data set in 2009 consists of two observations on Aug 24*^ and Aug 25*^ that 
are combined to increase the final [u,v] coverage. Given a rotation period of 55.15 days, 
starspots will migrate across the surface by ~6° over one night, so the combination of these 
two nights is not believed to adversely affect the quality of the extracted properties. Each 
night is the combination of observations using both the S1-E1-W1-W2 and S2-E2-W1-W2 


telescope arrays. Fig A-4 shows the distribution of the 704 [u,v] points obtained for the 
merged pair of observations. 

Season 2009 spans 130.7 days or ~2.4 rotation periods and has a AV = 0.154 mag. 
From one rotation to the next, the starspot properties do not appear to change significantly 
as illustrated by the low scatter compared to observation errors in Fig|^ 

Fig A-5 clearly show non-zero closure phases at both the lower and higher sampled 
spatial scales. An unspotted model image does not fit well with the measured interferometric 
data resulting in a reduced = 5.9. 


The best-fit model (reduced x^ = 1.44) contains three cool starspots. Fig A-5 contains 
the best-fit model image along with the model fits to the visibilities, triple amplitudes and 
closure phases. The starspot properties are listed in Table Al-1[ Fig |A-6| contains the final 
model, reconstructed and simulated images for the 2009 data set. The modeled starspot on 
the eastern limb was not conclusively detected in the reconstructed image. The properties of 
the two reconstructed starspots are nearly identical to the corresponding modeled starspots, 
to within errors. The western modeled and reconstructed starspots are included in the final 
results despite the measured covering factor of both being close to or below the CHARA 
array’s angular resolution (0.4 mas or 0 = 2.1%). The potential starspots are accepted as 
the model reduced is worse without its inclusion and the reconstructed starspot ads is 
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3.64. 


Two bright starspots are visible on the western limb; each starspot lies nearly eqnidistant 
above and below a east-west “equator”. These starspots are not visible in the simulated 
image. This allows for the possibility these starspots are genuine surface features. 


Based on the phased time series (see Fig|^, the interferometric observations were taken 
near maximum brightness. However as discussed in § |A-1 it is possible for starspots to be 
present on the visible surface even at maximum brightness. Given a rotation east to west, 
the identihcation of a large, cool starspot on the eastern limb would cause a drop in observed 
flux, as indicated by the photometry, as it rotates into view. 


The observing strategy of combining two consecutive nights of data provides a near 
500% gain in the best [u,v] coverage obtained in 2008. This has produced a much improved 
consistency between the model image and the reconstructed image along with a better quality 
of £t in both cases. However, the lack of multiple epochs does not provide consistency checks 
or a measure of the rotation period. 






1.38 0 - 1.38 


u (*10® rad'^) 

Fig. A-1.— The [u,v] coverage obtained for the 2008 observing run. A - Aug 17*^; 2010: B 
- Aug 18*^; C - Aug 19*^; D - Aug 20*^; E - Aug 21^*; F - Sep 20*^; G - Sep 27*^ 
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Spatial Frequency (rad'^) 

Fig. A-2.— The observed closure phases for the Aug 2008 data sets. Red Cross: Aug 
Ijth^ Asterisks: Aug 18*^. Yellow Squares: Aug 19*^. Green Diamonds: Aug 20*^. 

Blue Triangles: Aug 21^h The distinct non-zero closure phase signature points to surface 
asymmetries. The similarity in the closure phase between nights indicates a consistent 
asymmetric surface pattern from night to night. 
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Fig. A-3.— Results from the Aug 17*^, Aug 18*^, and Aug 19*^, 2008 data sets, including the 
model images (top row), reconstructed images (middle row), and simulated images (bottom 
row). The white dot in the lower right hand corner represents the 0.4 mas resolution limit 
for the CHARA array. 
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Fig. A-5.— The best fit results for the Aug 24*^ + Aug 25^^, 2009 data sets. Top Left: The 
model image. Top Right: The observed minus modeled visibilities as a function of baseline. 
Bottom Left: The closure phase as a function of spatial frequency. The orange asterisks 
indicate observed data and the red diamonds are the modeled £t. Bottom Right: The triple 
amplitudes as a function of spatial frequency. The symbols mean the same as in the closure 
phase plot. 
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Fig. A-6.— Results from the Aug 24*^ + Aug 25*^, 2009 data set, including the model images 
(left), reconstructed images (middle), and simulated images (right). The white dot in the 
lower right hand corner represents the 0.4 mas resolution limit for the CHARA array. 
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Table 1. 

CHARA Observing Log 

Date** 

Baselines 

# of [Ujv] points 

Calibrators*^ 

2008 

Aug 17*'* (54695.5) 

S1-E1-W1-W2 

96 

37 And (2), 45 Per (3) 

Aug 18*'* (54696.5) 

S1-E1-W1-W2 

144 

7 Lyr, 7 And (2), 37 And, C Per (2) 

Aug 19*'* (54697.5) 

S1-E1-W1-W2 

48 

7 And, C Per (2) 

Aug 20*'* (54698.5) 

S1-E1-W1-W2 

96 

7 And (2), 37 And (2), 45 Per (3) 

Aug 21«* (54699.5) 

S1-E1-W1-W2 

96 

7 And (2), 37 And (2), C Per, 45 Per 

Sep 20*** (54729.5) 

S1-E1-W1-W2 

72 

7 And (2), C Gas, 5 Per (2) 

Sep 27*** (54736.5) 

S1-E1-W1-W2 

72 

a Cyg, 37 And (2), ^ Per (2), tet Gem (3) 

2009 

Aug 24*'* (55067.5) 

51- E1-W1-W2 

52- E2-W1-W2 

272 

7 And (3), 37 And (2) 

Aug 25*'* (55068.5) 

51- E1-W1-W2 

52- E2-W1-W2 

432 

7 And (4), 37 And (2), HR 75 

2010 

Aug 2"** (55410.5) 

51- E1-W1-W2 

52- E2-W1-W2 

168 

7 And (2), 37 And 

Aug 3*''* (55411.5) 

51- E1-W1-W2 

52- E2-W1-W2 

456 

a Cyg, 7 And (3), 37 And (2) 

Aug 10*'* (55418.5) 

51- E1-W1-W2 

52- E2-W1-W2 

432 

cr Cyg, 7 And (3), 37 And (4) 

Aug 11*'* (55419.5) 

51- E1-W1-W2 

52- E2-W1-W2 

288 

cr Cyg, 7 And (4), 37 And (2) 

Aug 18*'* (55426.5) 

51- E1-W1-W2 

52- E2-W1-W2 

432 

a Cyg, 7 And (3), 37 And (5) 

Aug 19*'* (55427.5) 

51- E1-W1-W2 

52- E2-W1-W2 

432 

cr Cyg (2), 7 And (5), 37 And (7) 

Aug 24*'* (55432.5) 

51- E1-W1-W2 

52- E2-W1-W2 

528 

cr Cyg (2), 7 And (6), 37 And (6) 

Aug 25*'* (55433.5) 

51- E1-W1-W2 

52- E2-W1-W2 

384 

cr Cyg (2), 7 And (5), 37 And (2) 

Sep 2"** (55441.5) 

51- E1-W1-W2 

52- E2-W1-W2 

528 

7 And (7), 37 And (6) 

Sep 3*'*' (55442.5) 

51- E1-W1-W2 

52- E2-W1-W2 

600 

7 And (9), 37 And (3) 
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Table 1—Continued 


Date'^ 

Baselines 

# of [u,v] points 

Calibrators'^ 

Sep 10^'* (55449.5) 

51- E1-W1-W2 

52- E2-W1-W2 

336 

7 And (6), 37 And (2) 


2011 


Sep 2"^ (55806.5) 

S1-S2-E1-E2-W1-W2 

360 

O' Cyg, 7 And (2), 22 And (3), HR 653 

Sep 6‘'* (55810.5) 

S1-S2-E1-E2-W1-W2 

392 

a Cyg, 7 And (2), 22 And (3), HR 653 

Sep 10*'* (55814.5) 

S1-S2-E1-E2-W1-W2 

360 

7 And (2), 22 And 

Sep 14*'* (55818.5) 

S1-S2-E1-E2-W1-W2 

864 

7 And (4), 22 And 

Sep 19"* (55823.5) 

S1-S2-E1-E2-W1-W2 

808 

7 And (3), 22 And, HR 653 (2) 

Sep 24*'* (55828.5) 

S1-S2-E1-E2-W1-W2 

200 

7 And, 22 And, HR 653 (2), rj Aur 


^The number in parenthesis is the approximate MJD 

'^The number in parenthesis is the number of observations during the night 


Table 2. Calibrator Angular Diameters 


Calibrator (mas) 


Reference 


37 And 

0.682±0.03 

45 Per 

0.41±0.02 

7 Lyr 

0.740±0.100 

7 And 

0.663±0.024 

C Per 

0.703±0.021 

^ Cas 

0.280±0.020 

S Per 

0.555±0.038 

O' Cyg 

0.542±0.021 

tet Gem 

0.820±0.030 

HR 75 

1.04±0.012 

22 And 

0.591±0.041 

HR 653 

0.646±0.045 

r) Aur 

0.419±0.063 


new CHARA/MIRC measurement 
I Barnes et al.|l978 \ 
new CHARA/MIRC measurement 
new CHARA/MIRC measurement 
new CHARA/MIRC measurement 
new CHARA/MIRC measurement 
new CHARA/MIRC measurement 
new CHARA/MIRC measurement 
new CHARA/MIRC measurement 


CHARM2 

SearchCal i 

Richichi et al. 

20051 

Bonneau et al. 

2006 

SearchCal i 

Bonneau et al. 

2006 


new CHARA/MIRC measurement 
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Table 3. Photometry Observing Log 


Season 

Dates 

(MJD) 

To 

(MJD) 

AT 

(days) 


Period 

(days) 

AV 

(mag) 

2008 

54729.7979 - 54861.6223 

54722.0 

131.8 

66 

54.27±0.032 

0.165 

2009 

55091.9383 - 55222.6037 

55098.0 

130.7 

58 

55.15±0.91 

0.154 

2010 

55459.7624 - 55581.6084 

55488.7 

121.8 

73 

53.35±1.1 

0.099 

2011 

55823.6466 - 55958.5967 

55813.0 

134.9 

194 

53.30±1.9 

0.057 
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Table 4. 2010 Starspot Properties 


MODEL 

Param. 

Epoch 1 

Epoch 2 

Epoch 3 

Epoch 4 

Epoch 5 

Epoch 6 

<Pi (%) 

5.3±5.3 

7.6±4.0 

8.4±5.0 

7.9±3.7 

13.8±4.6 

44±12 

h (°) 

-10.1±1.3 

0.5±1.3 

7.1±1.2 

11.5±1.1 

4.9±1.017 

55±0.98 

h (°) 

-55.12±0.84 

-59.4±1.0 

-64.9±5.4 

-27.17±0.79 

-46.8±1.8 

-40.9±5.6 

Tri 

0.926±0.013 

0.758±0.018 

0.768±0.016 

0.772±0.036 

0.870±0.059 

0.925±0.017 

Starspot 

C 

E 

F 

F 

G 

— 

<f>2 (%) 

20.5±3.8 

5.8±3.1 

7.4±5.1 

7.1±4.6 

8.0±8.0 

13.5±6.4 

b2 (°) 

23.5±1.3 

-l.litl.l 

16.4±1.1 

34.9±1.4 

31.4±1.4 

29.84±0.96 

h (°) 

3.12±0.75 

-19.0±1.6 

-16.4±1.7 

19.00±0.96 

24.0±1.4 

9.7±4.9 

Tr2 

0.893±0.007 

0.860±0.015 

0.758±0.019 

0.756±0.057 

0.790±0.059 

0.898±0.048 

Starspot 

B 

D 

E 

E 

F 

G 

03 (%) 

11.5±5.5 

5.0±2.6 

5.3±3.4 

6.3±4.4 



h (°) 

52.5±2.1 

-6.8±1.0 

ll.Oitl.O 

28.3±1.5 

— 

— 

h (°) 

76.4±6.9 

10.6±1.2 

30.4±1.7 

70.8±7.8 

— 

— 

Tr3 

0.801±0.042 

0.860±0.015 

0.853±0.014 

0.853±0.0105 

— 

— 

Starspot 

A 

C 

D 

D 

— 

— 

04 (%) 


21.8±5.8 

4.0±4.1 




h (°) 

— 

54.95±0.99 

11.8±2.8 

— 

— 

— 

k (°) 

— 

77.5±1.5 

68.1±6.7 

— 

— 

— 

Tra 

— 

0.907±0.011 

0.853±0.017 

— 

— 

— 

Starspot 

— 

B 

C 

— 

— 

— 

"^red 

4.92 

2.87 

2.88 

4.14 

5.35 

4.79 

SQUEEZE 

01 (%) 




4.8 (-3.1) 

5.8 (-8.0) 


fel (°) 

— 

— 

— 

3.4 (-8.1) 

4.6 (-0.3) 

— 

h (°) 

— 

— 

— 

-28.74 (-1.57) 

-34.2 (12.6) 

— 

Tri 

— 

— 

— 

0.724 (-0.048) 

0.817 (-0.053) 

— 

^ds 

— 

— 

— 

11.55 

8.35 

— 

Starspot 

— 

— 

— 

F 

G 

— 

02 (%) 

9.0 (-11.5) 

4.0 (-1.8) 

5.76 (-1.64) 

4.8 (-2.3) 

5.8 (-2.2) 

6.8 (-6.7) 

h (°) 

23.6 (-0.1) 

-2.3 (-1.2) 

16.3 (-0.1) 

28.7 (-6.2) 

27.4 (-4.0) 

21.10 (-8.74) 

l2 (°) 

2.5 (-0.62) 

-27.4 (-8.4) 

-22.0 (-5.6) 

7.86 (-11.14) 

11.7 (-12.3) 

0.0 (-9.7) 

Tr2 

0.856 (-0.037) 

0.856 (-0.004) 

0.921 (0.163) 

0.724 (-0.032) 

0.761 (-0.029) 

0.849 PO.049) 

^ds 

4.12 

7.11 

16.20 

12.17 

9.81 

5.46 

Starspot 

B 

D 

E 

E 

F 

G 

03 (%) 


4.0 (-1.0) 

3.6 (-1.7) 




h (°) 

— 

-11.5 (-4.7) 

4.6 (-6.4) 

— 

— 

— 
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Table 4—Continued 


MODEL 


Param. 

Epoch 1 

Epoch 2 

Epoch 3 

Epoch 4 

Epoch 5 

Epoch 6 

fe n 

— 

9.4 (-1.2) 

28.8 (-1.6) 

— 

— 

— 

Tr3 

— 

0.856 (-0.004) 

0.849 (-0.004) 

— 

— 

— 

'^ds 

— 

6.28 

8.90 

— 

— 

— 

Starspot 

— 

C 

D 

— 

— 

— 

4>4 (%) 



2.6 (-1.4) 




b4 (°) 

— 

— 

-30.0 (-41.8) 

— 

— 

— 

i4 n 

— 

— 

2.7 (-65.4) 

— 

— 

— 

Tra 

— 

— 

0.915 (0.062) 

— 

— 

— 

'^ds 

— 

— 

6.16 

— 

— 

— 

Starspot 

— 

— 

C 

— 

— 

— 

Xled 

1.00 

1.00 

1.03 

0.98 

0.97 

0.99 


Note. — The parenthetical number beside the SQUEEZE parameter values is the differ¬ 
ence between the SQUEEZE and model values for the same starspot in each image. For 02 
and Tr values, a positive number indicates a larger value for the SQUEEZE parameter. For 
b and a positive number indicates the starspot in the SQUEEZE image is further north 
and west, respectively. 
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Table 5. Evidence of Stellar Rotation in the 2010 Data Set 


Starspot 

Epoch Range 

A0 

ATfl 

P rot 





(days) 

A 

1^.? 

— 

— 

— 

B 

1-5>2 

1 

-0.014 

46.6 

C 

1-5>2 

-0.3 

-0.066 

44.3 


2^>3 

-1 

-0.007 

47.9 

D 

2^>3 

-0.5 

-0.007 

56.9 


3^-4 

1 

0.000 

70.1 

E 

2^>3 

-0.2 

0.000 

63.7 


3^-4 

-0.3 

-0.002 

78.7 

F 

3^>4 

-0.5 

0.014 

77.0 


4-5>5 

0.1 

0.018 

63.4 

G 

5-5>6 

0 

0.028 

49.1 



Table 6. 2011 Starspot Properties 


MODEL 

Param. 

Epoch 1 

Epoch 2 

Epoch 3 

Epoch 4 

Epoch 5 

Epoch 6 

01 (%) 

16.9±5.7 

12.5±4.8 

11.8±3.1 

10.2±2.5 

10.4±6.4 

14.5±3.8 

fei (°) 

-1.3±2.1 

4.15±0.93 

11.9±1.3 

24.9±0.6 

1.30±0.92 

14.6±1.6 

h n 

-60.3±5.4 

-35.6±1.6 

-7.4±1.0 

11.8±1.3 

-60.7±3.5 

-34.1±1.2 

Tri 

0.825±0.014 

0.825±0.015 

0.797±0.013 

0.860±0.007 

0.825±0.011 

0.825±0.015 

Starspot 

B 

B 

B 

B 

C 

C 

02 (%) 

10.0±4.5 

10.0±5.4 

14.7±5.7 


14.8±5.3 

15.6±6.0 

b2 (°) 

3.65±0.95 

20.1±1.7 

37.7±1.4 

— 

40.5±1.4 

62.1±2.6 

fa n 

8.58±0.76 

38.1±1.2 

68.1±2.2 

— 

49.8±2.9 

85.6±5.7 

Tr2 

0.864±0.013 

0.864±0.015 

0.867±0.015 

— 

0.842±0.014 

0.849±0.016 

Starspot 

A 

A 

A 

— 

B 

B 

"^red 

2.91 

3.29 

3.81 

4.13 

6.24 

5.69 

SQUEEZE 

01 (%) 


4.0 (-8.5) 

7.8 (-4.0) 

4.8 (-5.4) 

4.0 (-6.4) 


bi n 

— 

11.54 (7.39) 

16.3 (4.4) 

24.8 (-0.1) 

-22.33 (-23.63) 

— 

h n 

— 

-29.3 (6.3) 

-14.5 (-7.1) 

6.3 (-5.5) 

14.4 (75.1) 

— 

Tri 

— 

0.839 (0.014) 

0.751 (-0.046) 

0.701 (-0.159) 

0.849 (0.024) 

— 

^ds 

— 

9.09 

5.59 

9.67 

3.23 

— 

Starspot 

— 

B 

B 

B 

C 

— 

02 (%) 

4.8 (-5.2) 



2.6 

4.8 (-10.0) 


b2 n 

9.21 (5.56) 

— 

— 

3.4 

24.8 (-15.7) 

— 

fa n 

12.88 (4.30) 

— 

— 

-31.4 

5.1 (-44.7) 

— 

Tr2 

0.860 (-0.004) 

— 

— 

0.915 

0.849 (0.007) 

— 

^ds 

5.31 

— 

— 

3.38 

3.18 

— 

Starspot 

A 

— 

— 

B 

B 

— 

"^red 

1.01 

1.02 

0.98 

0.99 

0.95 

0.97 


Note. — The parenthetical number beside the SQUEEZE parameter values is the difference between the 
SQUEEZE and model values for the same starspot in each image. For <p 2 ^.nd Tr values, a positive number 
indicates a larger value for the SQUEEZE parameter. For h and a positive number indicates the starspot in 
the SQUEEZE image is further north and west, respectively. 



Table 7. Evidence of Stellar Rotation in the 2011 Data Set 


Starspot 

Epoch Range 

A0 

ATfl 

P rot 



(YoJ 


(days) 

A 

1-5>2 

0 

0.000 

43.4 


2-5>3 

5 

0.003 

45.9 

B 

1-5>2 

-4 

0.000 

56.7 


2^>3 

-1 

-0.028 

49.8 


3^-4 

-2 

0.063 

64.6 


4^-5 

5 

-0.035 

51.1 


5^>6 

1 

0.007 

59.2 

C 

5-5>6 

5 

0.000 

61.2 


Table Al-1. 2009 Starspot Properties 


Aug 24*'"+25‘'" 


Parameter 


Model 

SQUEEZE 

01 

(%) 

16.0±1.9 

— 

bi 

n 

-8.6±2.0 

— 

h 

n 

-76.0±1.9 

— 

Tm 

(7 


0.961±0.008 

— 

02 

(%) 

4.1±6.8 

3.6 


n 

-2.3±1.9 

-3.4 

h 

n 

-13.4±5.3 

-26.8 

Tr2 


0.931±0.024 

0.922 

(7 


— 

4.8 

03 

(%) 

2.2±2.2 

2.0 

bs 

n 

-0.8±1.4 

-2.3 

h 

n 

22.9±2.0 

21.1 

Tr3 


0.979±0.010 

0.924 

(7 


— 

3.6 

Reduced 


1.44 

0.93 



